Transport through open quantum dots: making semiclassics quantitative 
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We investigate electron transport through clean open quantum dots ("quantum billiards"). We 
present a semiclassical theory that allows to accurately reproduce quantum transport calculations. 
Quantitative agreement is reached for individual energy and magnetic field dependent elements of 
the scattering matrix. Two key ingredients are essential: (i) inclusion of pseudo-paths which have 
the topology of linked classical paths resulting from diffraction in addition to classical paths and 
(ii) a high-level approximation to diffractive scattering. Within this framework of the pseudo-path 
semiclassical approximation (PSCA), typical shortcomings of semiclassical theories such as violation 
of the anti-correlation between reflection and transmission and the overestimation of conductance 
fluctuations are overcome. Beyond its predictive capabilities the PSCA provides deeper insights into 
the quantum-to-classical crossover. 

PACS numbers: 



INTRODUCTION 



The ability to controllably fabricate, manipulate 
and examine structures on the sub-micrometer scale 
has let to the observation of quantum phenomena 
in electron transport such as, e.g., universal conduc- 
tance fluctuations (UCF) in chaotic billiards and weak 
localization (WL), which dominate transport at the 
nanoscale. • By reducing the characteristic system size 
below the electronic inelastic mean free path, transport 
enters the so-called ballistic regime. 3 Ballistic electron 
transport is a prime candidate for semiclassical descrip- 

y ong 4, 5, 6, 7, 8, 9, 10, 11, 12, 13, 14, 15, 16, 17, 18, 19, 20, 21, 22, 23, 24,25, 26 

where the classical trajectories carry an amplitude which 
reflects the stability of the classical orbits and a phase 
which contains the classical action and accounts for 
quantum interference. 2 On a more fundamental level, 
the semiclassical framework provides a conceptually 
powerful bridge between classical and quantum mechan- 
ics allowing an intuitive approach to quantum mechanics 
and quantum chaos in general, and to transport through 
open quantum dots or so-called quantum billiards in 
particular. ' 5 ' 6 

Several semiclassical approximations (SCAs) based 
on the approximation of the constant-energy Green's 
function for propagation in a billiard have been proposed 
and compared with numerical quantum transport cal- 
culations or experiment. 8 ^ 9 - 10 ' 11 ' 12 ' 13 - 14 ' 15 ' 16 ' 17 - 18 ' 19 ' 28 ' 29 
While many qualitative features could be well repro- 
duced, quantitative agreement on a system-specific level 
has remained a challenge. 

One underlying difficulty is the multi-scale nature of the 
quantum-to-classical transition for transport through 
open quantum dots. For the semiclassical approxi- 
mation to hold, the de Broglie wavelength A should 
be vanishingly small compared to all characteristic 



dimensions of the device. Such asymptotic theories 
have been successfully employed to reproduce, upon 
ensemble averaging, random matrix theory (RMT) re- 
sults for chaotic cavities (see e.g. Ref. 7,20,21,22,23). A 
quantitative comparison on a system-specific level with 
full quantum calculations or experiments is, however, 
only possible in the non-asymptotic regime where A is 
small compared to the linear dimension D of the dot, 
A <C D, but still comparable to the width of the lead 
(or quantum wire) d, A < d. Moreover, for billiards 
with sharp edges the proper asymptotic limit is, rigor- 
ously, out of reach. The present theory addresses this 
non-asymptotic semiclassical regime, often also referred 
to as the "near" semiclassical regime. We show that 
the proper inclusion of diffractive contributions allows 
to quantitatively reproduce quantum calculations. The 
diffractive coupling between classical paths gives rise 
to pseudo-paths that are missing in the standard SCA 
and are the key to remedy many of the deficiencies of 
semiclassical approximations. 

We show in the present communication that this pseudo- 
path semiclassical approximation (PSCA) can reach 
quantitative agreement with full quantum simulations 
provided a high-order diffraction theory for the coupling 
between classical paths is used. For the scattering at the 
leads we develop an approximation involving elements of 
both the uniform theory of diffraction (UTD) 30 - 31 and 
the geometric theory of diffraction (GTD) 32 referred 
to in the following as the GTD-UTD approximation. 
With these ingredients good agreement with quantum 
simulations is found. 

One key conceptual insight is the unambiguous identifi- 
cation of the paths that contribute to quantum transport. 
We apply the present theory to a circle-shaped regular 
quantum dot for which the enumeration of paths, more 
precisely of path bundles, is easily possible. Unlike for 
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FIG. 1: (Color online) Geometry of the circular billiard of 
radius p with perpendicular leads of equal width d. In ac- 
cordance with previous work 3,i ' 34, ' io on the circular and the 
stadium billiard of equal area a — 4 + n, we choose d — 0.25 
and p — yl + 4/ir. 

chaotic dots, for which the exponential proliferation of 
contributing paths as a function of the pathlcngth makes 
their unique identification difficult, for regular systems 
their enumeration and identification is straightforward 
up to large pathlengths. 

The circular billiard is depicted in Fig. 1. The leads 
are attached at right angle and have equal width d. 
In order to probe the local topology of the cavity we 
require sufficiently long dwell times such that d/p <C 1. 
The wavelength of the electron A fulfills A « p for the 
semiclassical limit to hold inside the cavity. However, 
as in experimental or numerical studies of quantum 
billiards A < d. Our semiclassical theory can thus 
be quantitatively compared with quantum mechanical 
numerical calculations for the circular billiard. 
This paper is organized as following: In Sec. II we 
review both the standard SCA as well as the PSCA. 
These approximations differ by the different path sets 
entering the corresponding Green's function. In Sec. Ill 
we introduce the GTD-UTD diffraction approximation 
which is a key to the quantitative agreement between the 
PSCA and quantum mechanics for transport properties. 
The important role of pseudo-paths and a proper diffrac- 
tion theory is demonstrated on the level of quantum 
mechanical length-area spectra 16 in Sec. IV. Finally, we 
compare in Sec. V the semiclassical predictions for a 
variety of quantum transport properties that play a key 
role in the understanding of the quantum-to-classical 
crossover, in particular conductance fluctuations (CF), 
weak localization (WL) and quantum (non-thermal) 
shot noise, with quantum calculations. 

II. QUANTUM TRANSPORT THROUGH 
BILLIARDS 

We consider ballistic transport through quantum bil- 
liards for which transport properties are determined by 
the wave number (k) and magnetic field (B) dependent 
quantum mechanical Hamiltonian S-matrix. Dissipative 



or dephasing processes are neglected. (We will return to 
the effect of decoherence below.) Moreover, we refer to 
a "clean" billiard when elastic scattering at a disorder 
potential in the interior of the structure is absent. In 
this limit, the S-matrix represents elastic scattering at 
the boundaries of the billiard only. 

The S matrix elements Sn'm (fc, B) describe the scatter- 
ing from transverse mode m in lead i to mode n in lead 
j. We denote the transmission amplitudes from lead I 
to lead 2 as t nr „(£;,i?) = Snm\k,B) and the reflection 
amplitudes back into lead I as r nm (k 1 B) = Snm\k, B). 
According to the Landauer formula the conductance g of 
a quantum billiard is directly proportional to the total 
transmission T(k, B), 

2 P 2 Op 2 JH, J£ 

g(k, B) = —T(k, S) = — EE I W*> B)\\ (1) 

n— 1 m— 1 

where M is the number of open modes in the leads having 
equal width d. The S matrix elements can be determined 
by a projection of the lead modes tp m (yi) (yi is the trans- 
verse coordinate in the lead) onto the Green's function 
G{yj,yi,k,B) for propagation from j/j to yj at the cav- 
ity lead junction of lead i and lead j, respectively. The 
S'-matrix elements at B = are given by the Fisher-Lee 
equations 36 

^im(^j B 0) i \J k X2 ^n^x\ 7 m X 

J dy 2 J dyi (p* n (y 2 ) G(y 2 , yi, k, B = 0) <p m (yi) , (2) 

where k Xl . m {k X2 . n ) is the longitudinal wave number in 
lead I (lead 2). The prefactor ^k X2tn k Xl>m is due to flux 
normalization. We use atomic units (ft = |e| = m e s = 1). 
At zero magnetic field the mode wavefunctions take the 
form 

, , [2 . miry 

y m {y) = y 2 s ~d~ ' ^' 

For non-zero magnetic field B ^ 0, Eqs. (2, 3) have to 
be modified (see, e.g., Ref. 33,34 and references therein). 
We use the modular recursive Green's function method 
to calculate the exact quantum mechanical S'-matrix el- 
ements for a given k and B (for details see Ref. 33,34). 

A. Semiclassical approximations to transport 

The starting point of semiclassical approximations are 
the S-matrix elements [Eq. (2)]. In a first step, the quan- 
tum mechanical Green's function G is replaced by a cor- 
responding semiclassical expression which represents the 
Fourier-Laplace transform in stationary phase approxi- 
mation of the semiclassical limit of the Feynman propa- 
gator. It describes propagation along classical paths with 
fixed energy. Depending on the class of paths included, a 
hierarchy of different semiclassical approximations to the 
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Green's function results. These are to be distinguished 
from the different level of additional approximations em- 
ployed in the evaluation of the integral [Eq. (2)] which 
projects the Green's function onto the lead function. The 
latter gives rise to another set of semiclassical approxi- 
mations to the S'-matrix elements. 

We focus first on the replacement of G by a semiclassical 
approximation. For the latter we consider two different 
levels of approximation, the standard semiclassical ap- 
proximations (SCA) and the pseudo-path semiclassical 
approximation (PSCA). Both result from the stationary 
phase approximation (SPA) to the full Feynman propaga- 
tor reducing the continuous set of paths entering Feyn- 
man's path integral to a discrete subset of paths. As- 
suming that the classical action S is much larger than h 
and that well localized and separated stationary points 
with 6 Si = exist, the standard SCA contains exclu- 
sively classical paths. However, near sharp edges of the 
cavity or near the cavity-lead junctions, the de Broglic 
wavelength is not negligibly small and the SPA will fail. 
This leads to diffractive corrections which can be taken 
into account within the framework of the pseudo-path 
semiclassical approximation (PSCA). One of its salient 
features is that the basic notion of a propagator consist- 
ing of a sum over a discrete set of paths is preserved. 
Diffraction effects simply appear as additional contribu- 
tions to the path sum [see Sec. II (c)]. 
The classical action for an electron moving along a path 
q is given by S q = kL q + Ba q /c, where L q is the length 
and a q is the directed enclosed area of the path. a q can 
have both positive and negative values depending on the 
rotational direction of the path. In all our semiclassi- 
cal calculations (standard SCA as well as the PSCA) the 
magnetic field enters only via the Aharonov-Bohm phase 
Ba q /c. The curvature of the paths as well as the effect 
of non-zero magnetic field on the diffraction at the lead 
(introduced in Sec. Ill) is neglected since we consider the 
regime of weak magnetic fields (p <C ck/B with p being 
the radius of the circular cavity). 



B. Standard semiclassical approximation 

The Green's function within the standard SCA enter- 
ing Eq. (2) is given by: 

2tt 



G SCA ( y] ,y t ,k,B) 



\l\ D i(y^Vi^ k )\ 



(27ri) ^ q . 

x exp[iSg(j/j,j/i,fc,B)-i^ 9 ],(4) 

where D q (yj,yi,k) is the deflection factor (a measure 
for the divergence of nearby trajectories) and fi q is the 
Maslov index. The deflection factor is defined as 



\Dg(yj,yi,k)\ = 



1 



d 2 s a 



dyjdy, 



(5) 



where Xi (yi) is the longitudinal (transverse) component 
of the trajectory's starting point (i) and end point (J), 
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FIG. 2: (Color online) Two-dimensional length vs. enclosed 
area distribution of classical paths within the open (black 
dots) and closed (orange dots) circular billiard for (a) reflec- 
tion and (b) transmission. Each point represents one classical 
path which connects the centers of each lead. The inset of (a) 
shows typical paths for the path topology of the three first 
branches. 



respectively. The Maslov index increases by two for ev- 
ery reflection at the hard wall boundary of the billiard 
and by one when passing a focal point along the trajec- 
tory. Eq. (4) contains a sum over all classical paths q 
connecting the entrance lead i with the exit lead j (see 
Sec. II Bl). 

Evaluation of Eq. (2) with G SCA , 



+ SCA 



(fc, -E>) — i\/kx2,nkx\,m X 



dyi ldy2V* n (y2)G SCA (y 1 , y2 ,k,B) l p m (y 1 ), (6) 



proceeds either numerically 11 or analytically by invoking 
another set of SPAs. It was recognized from very early 
on that the SPA as applied to Eq. (6) is poorly justified 
in the non-asymptotic regime when A < d. Therefore, 
various diffraction integral approximations have been 
proposed. 10,12 The transmission amplitudes then take the 
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form: 

x \J\D q (k)\exp [iS q (k,B) -i^fi q ] 

x Cm(6i,k,d), (7) 

expressed in terms of diffraction coefficients c m (0,k,d) 
describing the diffractivc coupling from the entrance lead 
mode m into the dot and from the dot into the exit lead 
mode n. Deviating from previous calculations, we intro- 
duce for c m ^ n (9, k, d) a combination of Keller's geometric 
theory of diffraction (GTD) 32 and the "uniform theory 
of diffraction" (UTD). 30,31 The derivation of c m (6,k,d) 
within the GTD-UTD is given in App. A 2. The inclusion 
of diffraction effects in terms of diffraction coefficients in 
Eq. (7) preserves the structure of the semiclassical trans- 
mission amplitude in terms of a discrete sum over paths 
contributing to transport. The diffraction coefficients 
provide 6- and fc-dependent weighting factors for each 
path contributing to the transmission from mode m to 
mode n. Within the framework of the diffractive cou- 
plings into and out of the leads, the entrance and exit 
leads are treated as point scatterers. 15 Each path bundle 
which connects the entrance and the exit lead is replaced 
by an appropriately weighted representative path q con- 
necting the center of the entrance lead with the center of 
the exit lead. Consequently one can replace the deflec- 
tion factor as given in Eq. (5) by its value in the closed 
circular billiard D q = 1/ ^JkL qi where L q is the length of 
the path. 

The diffractive lead-dot couplings in Eq. (7) should be 
distinguished from diffractive corrections included in the 
propagation in the interior of the billiard. We refer to 
Eq. (7) as the standard SCA while inclusion of diffrac- 
tive corrections in the billiard corresponds to the PSCA. 

1. Paths entering the standard semiclassical Green's 
function 

Classical paths in a regular billiard (such as the cir- 
cle) feature a highly ordered structure of their length 
and enclosed area distribution (see Fig. 2). The branch 
structure of the "length-area" distribution of paths con- 
necting the entrance with the exit point is a specialty 
of the circular cavity. Along each branch the number of 
bounces off the wall increases by one from one path to 
its next higher neighbor. The points of convergence of 
each branch mark those paths that bounce off the wall 
infinitely many times and thus run exactly along the cav- 
ity boundary. In the limit 9 — > tt/2 (where 6 is the en- 
trance angle as given in Fig. f ) each branch contains an 
infinite number of paths. In our numerical calculations 
(see Sec. V) paths near this cluster point effectively do 
not contribute as they are cut off by vanishing diffraction 
coefficients c n (6 — > tt/2, k) — > 0. 



The distribution of paths eventually reflected back to the 
entrance point [Fig. 2 (a)] is symmetrically distributed 
relative to the a = axis due to time-reversal symmetry. 
Every path has a counterpart of equal length but oppo- 
site sign of the enclosed area. The lowest left and right 
branches consist of polygons with a number of revolutions 
71r = I in the cavity. The polygons can be characterized 
by an angle <fi — ^zl where up is the number of corners. 
Along each branch nc increases by one from one path to 
the next. The paths of the next higher branches revolve 
twice (n K = 2) around the circle and 6 = ^lhjl = izL. 
All higher branches can be described analogously. 
The branches of transmitted paths are not symmetric rel- 
ative to the a = axis but show a clear off-set (Fig. 2 (b)). 
Path pairs with similar length do not have, in general, 
the same topology but differ in the number of bounces 
off the hard wall boundary. As a consequence, these path 
pairs have different Maslov indices and thus do not inter- 
fere constructively. Fig. 2 also illustrates the difference 
between the open and closed billiard, i.e., the effect of 
"path shadowing" or suppression of longer paths due to 
their prior exit from the structure. All paths that would 
be geometrically reflected off the closed lead in the closed 
system are missing in the open billiard. The difference 
between the closed and open billiard is particularly evi- 
dent in reflection since all paths with four-fold symmetry 
leave the cavity via lead 2 before being reflected back to 
lead I. This is system-specific for the circular billiard 
with perpendicular leads (see Fig. I). 
The distinctly different path distributions for transmis- 
sion and reflection point already to a clear structural de- 
ficiency of the SCA. Many quantum properties of trans- 
port are a consequence of the intrinsic coupling between 
transmission and reflection. The standard SCA does not 
incorporate this quantum aspect of non-locality. Classi- 
cally, the path sets of transmission and reflection are dis- 
joint. The distribution of transmitted paths is markedly 
different from the one of reflected paths. In quantum 
transport transmitted and reflected paths are intertwined 
and must share the information on the relative phases. 
A semiclassical theory that reproduces quantum features 
must therefore allow for coupling between the path sets 
associated with transmission and reflection. This is the 
key property of pseudo-paths discussed in the following. 



C. Pseudo-path semiclassical approximation 

The pseudo-path semiclassical approximation 14,15 goes 
beyond the standard SCA by systematically including 
diffractive corrections into the propagation in the inte- 
rior of the billiard. In the present case, diffractive cor- 
rections arise from multiple back scattering, i.e., internal 
reflections at the leads. We point to the conceptual sim- 
ilarity to "pseudo-orbits" introduced in Ref. 37 for the 
study of the density of states in a closed wedge billiard 
and to "Hikami boxes" 2,38 introduced for elastic scatter- 
ing at short-ranged potentials in the interior of a diffusive 
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quantum dot. 

In line with multiple scattering theory the pseudo-path 
scmiclassical Green's function can be expressed in terms 
of a (semiclassical) Dyson equation, 15 



gPSCA _ QSCA i QSCAyQPSCA 



g sca Y,(vg sca )\ 

i=0 



(8) 



In the present case, G SCA plays the role of the unper- 
turbed Green's function and the perturbation "poten- 
tial" V accounts for the (internal) diffractive scatterings 
at the lead opening. The unperturbed propagation inside 
the cavity G SCA is equivalent to the free propagation in 
the "closed" system. 

For a two-terminal system the perturbation potential V 
is given by 



V 



Vi 
v 2 



(9) 



where V\ and describe the diffractive scattering off 
lead 1 and lead 2, respectively The semiclassical Dyson 
equation, Eq. (8), reads 

r<PSCA ^>PSCA\ /r<SCA r i SCA\ 

"1,1 "1,2 \ _ I "1,1 "1,2 \ 

/-'PSCA ^>PSCA / ~~ \ r'SCA r<SCA J 
"2,1 "2,2 / \"2.1 "2,2 / 

/G?ca Gf% A \ f Vl 0\ /GH CA G^ A \ 
+ W 2 CA g sca) {0 V J ^jpsca G psca; ■ I > 

To first order, the PSCA to the Green's function (denoted 
by G PSCA ^) connecting lead i with lead j includes terms 
in V of the form. 



g psca(i) = G^ A V t G 



SCA 
Li 



1=1,2 



SCA 



(11) 



Eq. (f f) may serve as example to illustrate the physics 
entering Eq. (10). It describes propagation from lead i to 
lead j via one intermediate visit to lead I where diffractive 
internal back scattering with amplitude v(9 q > ,9 qn k,d) 
takes place. The path from lead i to lead I, qij as well 
as from lead / to lead j, q'j z , are classical paths described 
by G SCA . Diffractive scattering couples the incident path 
qi i with angle 9 qi to the exiting path q'^ l with angle 9^ 
thereby coupling two disjoint subsets of classical paths 
and generating a first-order pseudo-path. The determi- 
nation of the diffraction coefficient v{9 q ' i ,9 qn k,d) will be 
discussed in more detail in Sec. III. From a conceptual 
point of view, the pseudo-path semiclassical approxima- 
tion [Eq. (10)] is closely related to the diagrammatic per- 
turbation theory. 1 Both leads 1 and 2 (i.e., V\ and V 2 ) 
act much like "Hikami boxes", 2,38 where electrons can- 
not be described semiclassically since the characteristic 
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FIG. 3: (Color online) Length-area distribution of first-order 
pseudo-paths (blue dots) for reflection (a) and transmission 
(b). Zeroth-order pseudo-paths (i.e., the classical paths) are 
denoted by orange dots (same as Fig. 2). The number of first- 
order pseudo-paths up to length L — 40 is larger by more than 
a factor of 60 compared to the number of classical paths. 



potential length scale (in our case the sharp edges of the 
leads) is smaller than the electron wavelength A. Thus, 
the wave nature of electrons has to be taken into account 
as diffraction at the point scatterers allowing classically 
distinct path sets to mix. This is crucial for the quantum 
corrections to the transmission and reflection amplitudes 
and repairs some of the deficiencies of the standard SCA 
(see Sec. V). 



1. Paths entering the pseudo-path semiclassical Green's 
function 

Within the PSCA pseudo-paths are formed by joining 
classical paths together via diffraction. With each in- 
creasing order of the PSCA the length-area distribution 
gets more and more densely filled with paths, or cquiva- 
lcntly, the number of paths increases with the order of the 
PSCA (Fig. 3). Pseudo-paths form product sets of classi- 
cal paths, e.g., joining a given classical path for transmis- 
sion with a classical path contributing to reflection (or in 
reverse order) forms the set of transmitted pseudo-paths 
to first order. Reflected first-order pseudo-paths result 
from joining two classical transmitted paths or two clas- 



sical reflected paths. This coupling allows to recover the 
non-locality of quantum transport. Higher-order pseudo- 
paths are constructed analogously. With increasing or- 
der and increasing (combined) pathlcngth the total num- 
ber of pseudo-paths exponentially proliferates. This is in 
sharp contrast to the power-law growth of purely classical 
paths for regular systems and explains why the effect of 
diffractive scattering is more likely visible in regular than 
in chaotic systems where the exponential proliferation of 
classical orbits may mask the diffractive contributions. 
The length-area distribution of first order pseudo-paths 
contained in Eq. (10) contributing to reflection is (of 
course) still symmetric [Fig. 3 (a)]. In addition, new 
branches appear with classical and pseudo-path part- 
ners of approximately equal length. The change in the 
branch structure is more dramatic in the spectrum of 
transmitted paths [Fig. 3 (b)]. The pseudo-paths lead 
to a "symmetrization" of the length-area distribution. 
The symmetrization results primarily from paths which 
change their rotational direction through diffractive scat- 
tering [see the pseudo-path in the inset of Fig. 3 (b)]. 
Branches of classical paths arc now completed by sym- 
metric pseudo-path "partner" branches of approximately 
equal length and different enclosed area. 
The weight and the phase of interfering paths are strongly 
influenced by the diffraction coefficient v{9' ,6,k,d). As 
will be demonstrated in Sec. V, the previously employed 
Fraunhofcr theory of diffraction 12,13 ' 14,15 is not suffi- 
ciently accurate as to give quantitatively reliable results 
for transport properties. The same holds for the Kirch- 
hoff theory of diffraction 10 which is closely related to the 
Fraunhofcr theory of diffraction and gives similar results 
for the diffraction coefficients (see Fig. 5 in the follow- 
ing section). For a quantitative agreement between the 
PSCA and quantum mechanics, it is thus necessary to 
go beyond low-order diffraction approximations and use 
a more sophisticated theory of diffraction. We present 
such a theory in the following section. 

III. DIFFRACTION AT THE LEAD 

The contribution of a given classical path with path- 
length L q to the standard semiclassical Green's function, 
G SCA , is 

^fSCA gifeLg— i37r/4— i^fiq (12) 

q ~ ^2nkL q 

Eq. (12) is the basic building block entering the 
Dyson Eq. (10) together with the diffraction coefficient 
v(9', 8, fc, d). In line with the far-field approximation un- 
derlying diffraction theory, v(6' ,9,k,d) is assumed to be 
independent of the length of the path approaching or ex- 
iting the diffractive scattering region. 
A successful application of the pseudo-path semiclas- 
sical approximation outlined in the preceding subsec- 
tion requires the determination of accurate diffraction 
coefficients v(9' ,9,k,d). Different approximations have 
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FIG. 4: (Color online) Sketch of diffraction at an open lead 
mouth: a path q reaches the orifice under an angle 9 and is 
backscattered into a path q' that leaves with an angle 9' (the 
angles 9 and 9' as depicted in the figure have opposite signs). 
The dashed lines denote that r and r ' are in the far field 
region. 



been used in the past for the inclusion of diffraction 
effects: Kirchhof diffraction approximation (KDA), 10 
Fraunhofcr diffraction approximation (FDA), 12,13 ' 14,15 
geometric theory of diffraction (GTD), 32 and the "uni- 
form theory of diffraction" (UTD). 30,31 We have devel- 
oped a new theory for the reflection at open lead mouths 
by combining the GTD with the UTD (the GTD-UTD) 
to take into account paths that scatter multiple times 
between the edges of the leads, (see appendix A 1). With 
this theory we obtain diffraction coefficients in excellent 
agreement to quantum mechanics. 

Consider, as a test case, the diffractive scattering (Fig. 4) 
at the lead mouth described by the first-order term 
[Eq. (11)], 

G P q ^ A(1) =G s q P A v(9>Ak,d)G s q CA . (13) 

The incoming path q is incident at angle (measured 
with respect to the surface normal) and is diffractively 
scattered into angle 9' under which path q' leaves the 
scattering region. We compare (Fig. 5) the present GTD- 
UTD theory with the FDA, the KDA and exact quan- 
tum mechanical (QM) calculations. 10 Even for a typical k 
value in the low mode regime (k =2.5ir/d), the agreement 
between the GTD-UTD and the exact QM calculations 
is very good whereas both the FDA and the KDA display 
clear deviations from the QM values. While these devi- 
ations do not appear dramatic at first glance, they are, 
in fact, quantitatively very important as the diffraction 
coefficient enters the Dyson series [Eq. (10)] to all orders. 
Note, however, that the GTD-UTD would fail in the limit 
9 — > it/ 2 [as indicated by the kink at 9' = 9 = 3tt/8 in 
Fig. 5 (c)] . This deficiency is of no concern for the present 
applications as the probability for diffractive scattering 
tends to zero in this limit. The maximum of the diffrac- 
tive reflection probability is clearly around the specular 
value 8' ps —9. However, it is important to note that 
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FIG. 5: (Color online) Absolute square of the diffraction co- 
efficient \v(0' , 8, k, d)\ 2 for diffractive scattering (see Fig. 4) 
within the FDA, the KDA, the GTD-UTD and exact quan- 
tum mechanical data for angle of incidence (a) 8 = tt/8, (b) 
8 = 7r/4, and (c) 8 — 3n/8 at k — 2.5n/d. The vertical lines 
in each frame mark the angles of specular reflection. The 
quantum mechanical (QM) and KDA coefficients are taken 
from Ref. 10. Inset of Fig. (b): the real and imaginary part 
of v{6', 8, k,d) within the GTD-UTD. 



the probability distribution possesses a local maximum 
at the backscattering angle 9' ~ 6. A non-negligible part 
of the electron wave is back-scattered into the direction 
from where it came from. This diffractive back-reflection 
should not be confused with the well-known Andreev 
back-reflection in which the back-reflected particle simul- 
taneously undergoes a particle-hole exchange. 39 Back- 
reflected paths are responsible for the symmetrization of 
the distribution of transmitted paths [Fig. 3 (b)] and are 
crucial for the understanding of the weak-localization dip 
in the transmission. 16 

While the KDA and the FDA have been successfully used 
in the past to explain certain features of conductance 
fluctuations, 10 ' 12 ' 13 ' 14 ' 15 only the GTD-UTD is precise 
enough to reproduce transport semiclassically on a quan- 
titative level (see Sec. V). In Sec. VB we compare the 
results for transport properties obtained by implementing 
the GTD-UTD and the FDA to full quantum mechanical 
calculations. 




FIG. 6: (Color online) Two-dimensional distribution of 
the absolute square of the quantum S'-matrix element 
\t 2 2{k,B)\ 2 = T 22 (k,B) as a function of the wavenumber k 
and the magnetic field B. 



IV. PATHS IN QUANTUM TRANSPORT 

The information on paths governing quantum trans- 
port can be reliably extracted from the two-dimensional 
Fourier transforms of the quantum mechanical S-matrix 
elements, 16 S nm (k, B). The S-matrix elements display 
a strongly fluctuating pattern as a function of k and B 
(see Fig. 6 for T^iik^ B)). Since the canonically conjugate 
variables to the wavenumber k and the magnetic field B 
are the length L and the directed area a, respectively, 
the Fourier transform 

S^JL, a) = fdk fdB e-^ kL+ ^S^(k, B) (14) 



allows the unambiguous identification of quantum paths 
contributing to quantum transport via their length and 
enclosed area [Fig. 7 (a)-(d)]. No apriori-assumption 
as to the existence of classical paths q(L,a) with path- 
length L and area a enters Eq. (14). The two-dimensional 
pathlcngth-area spectrum allows to identify both classi- 
cal as well as non-classical contributions to the full quan- 
tum spectra. Fig. 7 displays examples of pathlength-area 
spectra |S y (L,a)| 2 . [Note that Fig. 7 (c) is the absolute 
square of the Fourier transform of the transmission am- 
plitude whose absolute square is plotted in Fig. 6]. Obvi- 
ously the strong fluctuations of conductance in quantum 
transport (Fig. 6) are the result of the interference of 
clearly identifiable (quantum) paths (Fig. 7). 
The quantum mechanical S'-matrix elements S nm (k,B) 
are determined with the help of the modular recur- 
sive Green's function method 33,3 and then numerically 
Fourier transformed. Finite discretized intervals must 
be used when performing the Fourier transform numer- 
ically. The integration intervals are denoted by Afc = 
k m ax - k rnin (AB = B max - B min ) and the numerical 
grid spacings by 6k {SB). Accordingly, the resolvable 
length is AL = 2n/5k and the resolvable area interval 
is Aa — 2ir/(8B/c). The magnitude of the S'-matrix 
elements decreases with increasing length, which is an 
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FIG. 7: (Color online) The absolute square of the Fourier transform |S(L,a)| 2 as a function of length L and enclosed directed 
area a: (a) i?22(A a )> (b) Rn(L,a), (c) T22(L,a), and (d) Tii(L, a). The color shading is determined by log (\S(L, a) 2 + 1). 
The insets show the probabilities along the entire resolved length. The integration of Eq. (14) is performed numerically over 
k £ [2.2, 3.45]-7r/d discretized with 251 points and B/c G [—3,3] with 121 points. The spectra are compared to the distribution 
of classical paths (black dots) at zero B field (the curvature of the trajectories in the present magnetic field and energy regime 
is negligible). 



obvious consequence of open systems: the probability to 
stay within the cavity decreases with increasing length. 
The parameter 8k must be chosen sufficiently small such 
that the maximum resolvable length (AL) lies already in 
the region of strongly reduced amplitudes. Otherwise, 
the Fourier spectrum is visibly back-folded onto the fun- 
damental interval. The magnetic field interval is further 
restricted by the requirement that the curvature of the 
paths is negligible, i.e., the cyclotron radius ck/B should 
be much larger than the circle radius p. We chose the in- 
terval Ak and AB as well as the number of interval points 
such that a maximum length of AL = 100 is resolved [in 
Fig. 7 (a)-(d) only the contributions with L < 40 are 
shown, the insets contain the entire spectrum]. Except 
for i?n, the absolute square of the S-matrix elements is 
considerably damped at a length of AL = 100. Thus, the 
graphs represent essentially the entire length-area spec- 
trum. Only for R\\, contributions with L > 100 are 
non-negligible which leads to a back-folding near L = 
[the lowest branch structure near L = 2 in Fig. 7 (b)] . 
The quantum length-area spectra provide detailed infor- 



mation on the paths contributing to transport through 
a specific system. They represent the paths entering the 
full Fcynman path integral. The following general trends 
can be observed: Long paths are more prevalent in S'- 
matrix elements connecting low mode numbers (in the 
present case they contribute more strongly to S'-matrix 
elements Sn than to £22)- Lower modes favor smaller 
entrance and exit angles that are associated with longer 
paths with a larger number of bounces off the cavity 
boundary. The most important observation is the re- 
markably close correspondence of the quantum mechani- 
cal length-area spectrum to its classical counterpart. Im- 
portant contributions are located near classical paths. 
Moreover the branch structure of classical paths is re- 
produced [Fig. 7 (a)-(d)]. On the other hand, there are 
distinct structures which do not correspond to classical 
paths and which can be identified using the distributions 
in Fig. 3 (a) and (b). The quantum mechanical length- 
area spectra confirm the existence and substantial role 
of non-classical paths, the pseudo-paths: these are those 
paths that are one or several times diffractively reflected 
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FIG. 8: (Color online) T22(L,a) as a function of length L and enclosed directed area a: (a) full quantum S'-matrix, (b) the 
PSCA using the GTD-UTD, and (c) the PSCA using the FDA for the diffraction coefficients. Several regions with non-classical 
paths are highlighted by red circles. 



off the lead before exiting the cavity. Two examples are 
given in the following: For T22 we identify an ensemble 
of diffractive paths which, among others, contribute to a 
symmctrization of the spectrum [non-classical branches 
in Fig. 7 (c), two of them together with the classical 
partner branches are marked by arrows]. i?n reveals the 
importance of paths that are geometrically reflected off 
the open lead (e.g., the periodic contributions near a = 
along the length axis belong to horizontal paths bounc- 
ing increasingly many times back and forth). 
The importance of a given class of paths to quantum 
transport can be delineated by inverting this decomposi- 
tion process. Deleting a selected class of paths (classical 
or non-classical) from Sij{L,a) and performing the in- 
verse Fourier transform gives rise to truncated S'-matrix 
elements Sij (k, B) from which certain path contributions 
have been removed in a controlled manner. This is the 
key to detailed quantitative tests of semiclassical the- 
ories. Since summation of the (P)SCA over arbitrar- 
ily long paths is prohibitively complicated we can com- 
pare truncated quantum and semiclassical S-matrix el- 
ements where both quantum and classical paths only 
up to a maximum pathlcngth L < L max are included. 
The length-area spectra also allow sensitive tests for the 
proper diffractive weight of a given class of paths in a 
semiclassical theory. To this end, we first calculate the 
S'-matrix elements within the PSCA and then perform 
the Fourier transform [Eq. (14)] in analogy to quantum 
calculations. To analyze the role of a proper diffraction 
coefficient we use either the GTD-UTD (which gives good 
agreement with quantum mechanics, see Fig. 5) or the 
FDA (with poor agreement with quantum mechanics, sec 
likewise Fig. 5). The fact that back-reflection into the 
cavity is poorly described within the FDA is mirrored in 
the semiclassical length-area spectrum where important 
non-classical (diffractive) contributions have a far too low 
weight [Fig. 8 (c)]. In particular, the diffractive change 
of the rotational direction is insufficiently described (see 
Fig. 5 for 8' > 0). A clear indication for the essen- 
tial role of the corresponding paths is the improvement 
within the GTD-UTD. The length-area spectrum within 



the GTD-UTD remarkably reproduces even fine details 
of the quantum mechanical spectrum [compare Fig. 8 (a) 
and (b)]. The non-classical path sets in Fig. 8 can be 
identified using the path distributions within the PSCA 
(Fig. 3). The length-area spectra do not leave any am- 
biguity as to which paths contribute by which phase and 
weight. By taking into account pseudo-paths and weight- 
ing them with the appropriate diffraction coefficients the 
quantum mechanical path spectrum is reproduced on a 
quantitative level. In the following section we demon- 
strate how the accurate representation of the length-area 
spectra within PSCA directly translates into quantitative 
reproduction of transport properties. 



V. APPLICATION TO TRANSPORT 
THROUGH REGULAR BILLIARDS 

Phase coherent ballistic transport governed by quan- 
tum interference influences the conductance in several 
important ways: the conductance strongly fluctuates 
as a function of the Fermi energy, the magnetic field 
or the cavity geometry (conductance fluctuations, CF). 
The conductance is, on average, suppressed compared 
to the classical prediction and increases with external 
magnetic field (weak localization, WL). The noise carries 
the signatures of quantum mechanical uncertainty (shot 
noise) . 

Path interference has been the key to the 
understanding of phase coherent ballistic 
transport ".8,9,10,11,12,13,14,15,16,17,18,19,20,21,22,23,24,25,26,40 

For chaotic systems and large mode num- 
bers quantum transport properties have been 
attributed 17 ' 18 ' 19 ' 20,21 ' 22,23 to the interference of clas- 
sically allowed paths (SCA). We demonstrate in the 
following for the circular billiard at low mode numbers 
(accessible to experiments) that CF, WL, and shot noise 
cannot be explained by the interference of classical paths 
alone. In this regime, the SCA overestimates the CF, 
does not reproduce the weak-localization dip and also 
shows poor quantitative agreement with the quantum 
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mechanical prediction for shot noise. These difficulties 
can be overcome with the PSCA when a high-level 
diffraction approximation, the GTD-UTD, is employed. 
For technical reasons, we perform in both the SCA and 
the PSCA the summation over paths only up to a maxi- 
mum cut-off length L max . Correspondingly, we truncate 
the full quantum scattering matrix elements by setting 
all elements S nm (L,a) with lengths exceeding L max to 
and carry out the inverse Fourier transform. This 
allows a quantitative comparison between semiclassical 
and quantum calculations unaffected by the (inevitable) 
unitarity deficiency of a truncated semiclassical path 
sum. 



A. Conductance fluctuations 

Conductance fluctuations (CF) have been identified as 
a direct manifestation of phase coherent transport (see 
e.g. Refs. 1,2,8,18,41,42,43,44 and references therein). 
The strong fluctuations of the conductance as a function 
of, e.g., energy originate from path interference and 
thus give evidence for the wave nature of electrons in 
quantum dots. The CF offer one of the most stringent 
testing grounds for a semiclassical theory, when good 
agreement on the level of each individual S'-matrix 
element is required. 

In the following we compare the results for CF within the 
PSCA with the GTD-UTD, the SCA and the quantum 
mechanical calculations as a function of k at vanishing 
magnetic fields B = 0. Results for B ^ averaged over 
k can be found in Sec. VB. 

The semiclassical and quantum mechanical results both 
display strong fluctuations of the conductance [i.e., the 
total transmission T(k,B = 0)]. Their amplitude is, 
however, extremely sensitive to any deficiencies in the 
semiclassical approximations (inappropriate weighting 
of paths, missing paths). We emphasize that the 
comparison of Fig. 9 is on a fully differential level. No 
energy or ensemble average is involved. Unsurprisingly, 
the agreement between the SCA and the quantum 
data is poor and on a level of qualitative agreement at 
best [Fig. 9 (a) and (b)]. The functional dependence 
of T(k, B = 0) and R(k, B = 0) seems only weakly 
related to the quantum mechanical prediction. By 
ensemble averaging (e.g., over a suitable k interval) 
these discrepancies would be (partially) removed (or 
masked). To correctly reproduce the fluctuations in 
T(k,B = 0) = EtlnKniKB = 0)| 2 , all individual 
mode-to-mode amplitudes t mn have to be accurate. 
Obviously, the contribution of pseudo-paths included 
in the PSCA but missing in the SCA significantly 
improves the agreement with the quantum conductance 
fluctuations [Fig. 9 (c) and (d)]. 

Even for averaged quantities, such as the averaged 
conductance (T) Ak and the variance of the conductance 
fluctuations a = \J (T 2 ) Ak — (T) 2 Ak , the deficiencies of 



the SCA are still visible, in particular at larger cut-off 
lengths Lmax- Both (T) Ak and o are overestimated 
(Fig. 10) as compared to quantum mechanics which can 
be attributed, in part, to the lack of correlation among 
purely classical paths. 

In contrast, the PSCA shows very good agreement with 
quantum mechanics for both the averaged conductance 
as well as for the variance. For averaged quantities 
(Fig. 10) the inclusion of up to 4 th order diffractive 
scattering for L max = 40 is sufficient while for fully 
differential quantities (Fig. 9) 5 th order corrections for 
L max = 40 still improve the agreement. 
We point out that for the present system, taking into 
account paths up to a length of L max = 50 which 
corresponds to 16 radial traversals through the billiard, 
one reaches 70% of the unitarity level within both the 
truncated quantum mechanics and the PSCA. The 
essential contributions to CF, WL and shotnoise (see 
the following sections) is thus rooted in this subset of 
relatively short paths. 



B. Weak localization 

Weak localization (WL) is a well-known quantum 
correction to the classical diagonal approximation (see 
e.g. Refs. 1,2,8,16,17,18,19,20,28,45,46,47,48,49 and ref- 
erences therein) where the later corresponds to the re- 
striction to terms q = q' in the double sum over paths 
when calculating |t^„ A | 2 from Eq. (7). Off-diagonal 
terms q q' give rise to quantum interferences which 
manifest themselves as an increase of the averaged total 
reflection (R(B)) Ak = Y, m ,n( R nm(B)) Ak at B = in 
form of a pronounced peak. Correspondingly, the aver- 
aged total transmission (T(B)) Ak = J2 m ,n( T nm(B)) Ak 
features a dip which is an immediate consequence of uni- 
tarity. For the investigation of the weak localization 
dip (peak) we employ an average over a small window 
Ak = [2.2 — 2.8]7r/d of the k dependence of the probabil- 
ities: (T nm (B)) Ak = J Ak dk\t nm (k,B)\ 2 , (R nm {B)) Ak = 
J Ak dk\r nm (k,B)\ 2 . 

Since semiclassical theories are, by construction, not 
necessarily unitary at a given level of approximation, 
the anti-correlated peak-dip structure near B = pro- 
vides a sensitive test for semiclassical approximations. It 
has been shown that the quantum anti-correlation be- 
tween reflection and transmission ((R(B)) Ak — (R(B = 
0))Afe) = -((T(B)) Ak - (T(B = 0)) A fe) requires a corre- 
lation of transmitted and reflected paths. 16 This correla- 
tion is absent in the SCA such that no transmission dip 
is reproduced [see Fig. 11 (b)]. 

The role of diffraction in a quantum billiard manifests 
itself by a very good agreement of the PSCA with quan- 
tum mechanical results [see Fig. 11 (a)]. Tests for WL 
as a function of the cut-off length L max of the paths as 
well as of the order of diffractive scattering included show 
that convergence toward the (truncated) quantum result 
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FIG. 9: (Color online) (Conductance) Fluctuations in the transmission T(k, B = 0) and reflection R(k, B — 0) as a function 
of fc. The maximal included length is L max = 40. (a) Reflection and (b) transmission within the SCA (red/green solid line), 
and quantum mechanics (black dashed line). For the diffraction coefficients c m (8,k,d) entering the SCA we have used the 
GTD-UTD. (c) Reflection and (d) transmission within the PSCA with the GTD-UTD (red/green solid line), and quantum 
mechanics (black dashed line). The PSCA includes diffractive scattering up to 5 th order. 




FIG. 10: (Color online) Comparison between the PSCA, the SCA, and quantum mechanics (QM) for (a) the averaged trans- 
mission (T)a/c as a function of the inverse cut-off length 1/L max and (b) the averaged variance of the conductance fluctuations 
a. The average is performed in the interval k G [2.2, 2.8]ir/d. The data within the PSCA are calculated in 3 rd order for lengths 
L max — 10 — 30, in 4 th order for L max = 40 and 5 th order for L max = 50. The QM results extend to the exact value (L max = oo). 
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FIG. 11: (Color online) Comparison between truncated quan- 
tum mechanics, the SCA, and the PSCA for the weak local- 
ization dip in transmission for different L ma x- (a) PSCA using 
the GTD-UTD for all diffraction coefficients (red solid line), 
(b) Standard SCA (red solid line) (using the GTD-UTD for 
the diffraction coefficients Cm(0,k,d) for entering and exiting 
the circular cavity.) Black dashed lines: quantum results. 



is reached for L max = 40 when diffraction up to fourth 
order is included. We note an improvement compared to 
previous third-order calculations 16 especially due to the 
inclusion of the (small) real part of the diffraction coeffi- 
cient v(6', 9, k, d). 

It is instructive to analyze the build-up of the weak local- 
ization peak in reflection and of the dip in transmission 
from individual S-matrix elements. In reflection only the 
diagonal elements Ru(B) and R22{B) show a peak while 
the off-diagonal elements Ri2(B) and R2\{B) exhibit a 
dip [Fig. 12 (c), (d)]. This is due to the fact that time- 
reversal symmetric paths contributing to R nm (B) with 
different parity of modes m and n acquire an additional 
phase-shift of tt. As an example consider R2\(B), enter- 
ing in mode m = 1 and exiting in mode m = 2 produces a 
phase shift of tt. Thus the time-reversal symmetric paths 
interfere destructively. 

In transmission the major contribution to a dip originates 
from T22{B) in the chosen energy window of two open 
modes [Fig. 12 (e), (f)]. Note that T 21 {B) = T 12 {B) be- 
cause of the Onsager relation tnm {k, B) = tmn (k, —B) 
and the symmetry of the circular billiard which ensures 

that tmn (k, — B) = tmn (k, B). Due to time-reversal 
symmetry, R(B) = R{—B) giving rise to a symmetric 
peak in R as a function of B (Fig. 12). For unitary trans- 
port, T(B) = 1 — R(B), which implies a symmetric dip 
in transmission. For the truncated quantum mechanics 
and semiclassics where long paths are omitted unitarity 
is not preserved and T(B) is not exactly symmetric but 
features a slight shift of the minimum toward B < 0. 
Overall, the agreement between the PSCA and the full 
quantum calculation (truncated at the same pathlcngth) 
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FIG. 12: (Color online) Weak localization for I/ ma x = 40 as 
produced by quantum mechanics: left column (a),(c),(e), and 
PSCA with the GTD-UTD: right column (b), (d), (f). 



is remarkable (Fig. 12). The residual minor deviations 
are mainly ascribed to deficiencies in the diffraction co- 
efficients for which we use an analytical far-field approxi- 
mation (Sec. Ill and appendix Al and A 2). As, e.g., the 
diffraction coefficients c m {9,k,d) enter R nm , T nm to the 
fourth power, a small deficiency in c m (9,k,d) can have 
a sizeable effect on i? nm , T nm . The offset of the total 
reflection R and transmission T within PSCA [Fig. 12 
(b)] compared to the quantum mechanical result [Fig. 12 
(a)] mainly originates from imperfections in R22 and T22 
[compare Fig. 12 (c), (d), and (e), (f)]. The imperfec- 
tion could be cured by including a correction factor of 
w 0.97-0.99 in c 2 (9 7 k,d). 

To demonstrate how sensitively transport properties de- 
pend on the weights of the contributing paths we show 
the results for weak localization calculated within PSCA 
but now using the FDA instead of the GTD-UTD for the 
diffraction coefficients (Fig. 13). (We have used 3 rd or- 
der of the PSCA here, since the deviation between 3 rd 
and 4 th order is much smaller than the errors due to the 
simpler diffraction theory.) It is striking that especially 
both R22 and T22 seem to be underestimated which is 
due to deficiencies in C2(9,k,d) within the FDA. It is 
worthwhile pointing out that the FDA diffraction could 
not be "repaired" by a correction factor, as the errors 
in R22 and T22 are different. In other words, transmis- 
sion is more "sensitive" to a correct implementation of 
pseudo-paths than reflection. Of particular importance 
are paths which change the rotational direction. 16 Since 
the FDA does not give sufficiently weight to this class of 
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FIG. 13: (Color online) Weak localization for I/ max = 40. 
Left column (a),(c),(e): quantum calculations, right column 
(b),(d),(f): PSCA with diffraction coefficients from FDA. 
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(see e.g. Refs. 1,18,21,42,50,51,52,53,54,55,56,57 and ref- 
erences therein). At zero temperature (T = 0), the time- 
dependent current noise is due to the granularity of the 
electron charge and carries information about the wave 
vs. particle nature of charge transport. The Fano fac- 
tor F measures the amount by which the noise in phase 
coherent transport is suppressed relative to the Poisso- 
nian value of uncorrelatcd classical electrons. Within the 
Landaucr-Buttikcr picture, F can be expressed as 50 



F 



(Tr{tHr ] r)) Ak _ (J2 n T nVn) Ak 



(Tr(tH)) 



Ak 



(Sri T n)Ak 



(15) 



with r n , r\ n being the eigenvalues of the Hermitian ma- 
trices t't and r^r, respectively. 

Calculating the shot noise power from non-unitary scat- 
tering matrices is obviously a delicate matter, as replace- 
ment of r*r by 1 — tH leads, unlike for unitary descrip- 
tions, to different results. Furthermore, such a replace- 
ment may result in negative and thus unphysical values 
for the shot noise power, as non-unitary scattering ma- 
trices allow for the possibility of having t„ > 1 such that 
(1 — t„) < (for very high mode numbers as, e.g., in 
Ref. 18 such a situation may, however, be unlikely). By 
using Eq. (15), such difficulties can be avoided as both 
the transmission and reflection eigenvalues r„, ?/„ are, by 
construction, real and positive not only for the truncated 
quantum calculation but also for the PSCA and the stan- 
dard SCA. The standard SCA result for F strongly devi- 
ates from the quantum mechanical data (see Fig. 14). For 
small cut-off lengths L max , the value of F is smaller but 
increases more rapidly than the quantum mechanical re- 
sult with L max . The PSCA yields very good agreement 
with the quantum mechanical result for the shot noise 
Fano factor F. Note that for the largest cut-off length 
L mEUi = 50 the Fano factor F is already converged to its 
asymptotic value F w 0.28 suggesting that long paths do 
not play a significant role for F . This result agrees with 
the finding 51,52 that the shot noise power is of similar 
magnitude for regular and chaotic billiards as differences 
in the dynamics are most strongly felt by very long paths. 



1/L„ 



FIG. 14: (Color online) The fc-averaged Fano factor F as a 
function of the inverse cut-off length Lmax. The average is 
performed over the interval k G [2.2, 2.8]ir/d. Comparison 
between PSCA with the GTD-UTD, the SCA, and quantum 
mechanics (QM). (L max = oo corresponds to the exact result.) 



pseudo-paths (see Fig. 5) the transmission dip cannot be 
well reproduced. 



C. Shot noise 

Another quantity characteristic for quantum trans- 
port is the quantum shot noise power of the current 



VI. SUMMARY 

We have presented a semiclassical theory which is able 
to quantitatively reproduce full quantum results for scat- 
tering through microstructures with a specific geometry, 
in the present case a circular shaped billiard with leads 
oriented 90° degrees relative to each other. The present 
approach does not invoke the limit of large mode num- 
bers, (where the de Broglie wavelength A is small rela- 
tive to the lead width) but requires A to be small only 
on the scale of the linear dimension of the microstruc- 
ture (the circle). This non-asymptotic semiclassical the- 
ory allows a direct comparison with quantum calcula- 
tions as well with experiments on a system-specific level 
for individual S'-matrix elements avoiding any ensemble 



14 



averaging or fit parameters. This level of agreement al- 
lows us to perform detailed semiclassical investigations of 
quantum transport quantities such as the conductance 
fluctuations, the weak localization, and the shot noise. 
Our studies show unambiguously that for reproducing 
these quantities correctly, two major ingredients are, in- 
deed, crucial: (1) the inclusion of "pseudo-paths" in the 
semiclassical propagator which arc diffractively backscat- 
tered from the interior side of the cavity openings, (2) a 
sufficiently accurate description of the diffraction coef- 
ficients for the injection, ejection and back-reflection of 
particle flux at the cavity openings. We meet the lat- 
ter requirement by developing a combined geometric and 
uniform theory of diffraction (GTD-UTD). Pseudo-paths 
are crucial for reproducing the conductance fluctuations 
in transport and lead to a reduction of its variance. Also 
for the weak-localization effect we find that pseudo-paths 
are crucial, as no signature of weak localization appears 
in the transmission through the circular cavity without 
their contribution (even when the advanced diffraction 
theory is employed for all truly classical paths). For the 
shot noise power we showed that a standard semiclassical 
calculation (without pseudo-paths) gives sizeable discrep- 
ancies. The inclusion of pseudo-paths leads to agreement 
with the quantum mechanical result. We emphasize that 
the parameter regime in which we have identified the 
above effects of pseudo-paths coincides with the typical 
situation in quantum transport experiments. The latter 
usually feature only a few open lead modes M. 
The present results raise several interesting questions for 
future investigations: the comparison between the semi- 
classical approximations (PSCA and SCA) and full quan- 
tum calculations were performed for truncated path sums 
up to a finite path length L < L max . The primary rea- 
son for the truncation was technical, as the number of 
diffractive pseudo-paths exponentially proliferates with 
L — > oo also for classically regular structures and ex- 
act path sums become prohibitively difficult to perform. 
There is, however, a second conceptual motivation. In 
the experiment, dccohcrcncc due to inelastic scattering 
limits phase-coherent transport to pathlcngths L < 1$, 
where the phase-decoherence mean free path 1$ typically 
allows only a moderate number of traversals across the 
cavity. The latter restriction rules out that very long 
paths with L > contribute to quantum interference in 
the experiment, a feature which is naturally incorporated 
by way of the cut-off length L max in our semiclassical the- 
ory. Clearly, such long paths can still provide incoherent 
contributions. The present approach may thus contribute 
to a semiclassical unterstanding of decoherence effects in 
regular cavities. 28 

In the present treatment of diffractive scattering, both 
internal diffraction at the open lead mouth giving rise 
to pseudo-paths as well as the coupling between leads 
and cavity was performed for sharp edges. The weight 
of diffractive contributions can be changed by "round- 
ing off" the lead opening. An investigation of the de- 
pendence of the weak localization on the smoothness of 



the edges is currently underway. 64 The introduction of 
rounded corners has, however, another profound effect, 
apart from changing the weight of diffractive scatter- 
ing: an open circle with rounded edges of the leads is 
no longer regular but features a mixed phase space. This 
raises the question as to the interplay between diffractive 
scatterings at the lead opening and chaotic scattering in 
the interior of the billiard. For generic chaotic systems 
a number of alternative semiclassical theories has been 
proposed to explain universal features of quantum trans- 
port (see e.g. Refs. 17,18,19,20,21,22,23 and references 
therein). These theories typically employ an ensemble 
average and rely on a h — > limit which makes them 
complementary to the present system-specific approach 
for finite H. Bridging the gap between these two frame- 
works would be of great interest. One key ingredient 
would be to clarify the interplay between the diffraction- 
based pseudo-paths and the chaos-based correlated clas- 
sical path pairs (Richter-Sieber orbits 20 ). The relative 
weight for characteristic quantum transport effects car- 
ried by these two classes of paths when both are present 
(as in a chaotic billiard with sharp edges) remains an 
open question. 
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APPENDIX A: GEOMETRIC AND UNIFORM 
THEORY OF DIFFRACTION FOR THE 
COUPLING OF QUANTUM LEADS TO A 
BILLIARD CAVITY 

In a quantum billiard, the electron propagation in 
the leads is determined by a quantum mechanical wave, 
while the propagation inside the ballistic cavity can 
be described semiclassically, i.e., by propagation along 
classical trajectories with a quantum phase. The strong 
scattering effects that occur especially in the low mode 
regime at the orifices have been described in past works 
on transport through open billiards by the Kirchhoff 
diffraction approximation (KDA) 10 or on the level of the 
Fraunhofcr diffraction approximation (FDA). 12,13 ' 14 ' 15 
In the high-mode regime diffraction has been neglected 
altogether. 8 - 17 ' 18 . 19 . 20 - 21 - 22 - 23 . 24 - 26 Both the KDA and 
the FDA perform about equally well (see, e.g., Fig. 5). 
For the identification of pseudo-paths in the Fourier 
spectra of the conductance fluctuations, the KDA and 
the FDA have been sufficiently precise. However, in 
order to recover unitarity of the semiclassical S'-matrix 
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FIG. 15: (Color online) Propagation between two points f and 
r' in a semi-infinite plane in the presence of an open (lead). 
Contributions from the direct path (black dotted line), from 
specular reflection (blue solid line) , and from scattering at 
the orifice (red dashed line) are depicted. 



and weak-localization in transmission, a more precise 
diffraction theory must be implemented. For multiple- 
scattering paths, higher order products of the diffraction 
weights occur and even small errors are rapidly am- 
plified. We have therefore implemented a combination 
of the geometric theory of diffraction (GTD) 32 and 
the uniform theory of diffraction (UTD) 30 ' 31 for the 
diffraction coefficients in open billiards referred to in 
the following as the GTD-UTD. Both theories have 
been previously applied separately to the calculation of 
higher-order scattering corrections to Gutzwillcr's trace 
formula in closed quantum-billiards. 31 We discuss in 
section A 1 diffractive backscattering at a semi-infinite 
half-plane with an orifice and in section A 2 diffraction 
during propagation from a lead into a semi-infinite plane. 



1. The GTD-UTD for backscattering into the 
cavity 

The propagation between two points in a semi-infinite 
plane with a connected lead is depicted in Fig. 15. In the 
absence of a lead, the propagation between two points 
r and r ' in a semi-infinite plane is described by a sum 
of two Green's functions, G(f',f,k) = G dlr (r ', r, k) + 
G refl (r ', r, k). The first contribution corresponds to the 
direct path from r to r 1 (black dotted line) and is in the 
scmiclassical approximation [Eq. (12)] given by 



G dlr (f',r,fc) = 



^ikL(f\r)—i37r/4:—ifj l TT/2 



v / 2nkL(r',r) 
G SCA (r',f,fc), (Al) 



where L(r',r) is the distance between the two points 
and \i is the Maslov index (/.t = in this case). The 
second term corresponds to a propagation via a classical, 
specularly reflected path (blue solid line): 

ifc[L(r / ,a)+L(a,r)j— z3tt/4— iir 

G rcfl (r ', r, k) = -^^^^^^^^=. (A2) 
y/2nk [W.aj + L^r )] 



The orifice gives rise to an additional scattered wave for 
which we assume in the far-field limit (\kf\ 3> 1 and 
\kr'\ 1) a cylindrical wave emanating from the center 
of the orifice (point b) . Invoking far-field approximations 
is at the heart of semiclassical diffraction theories, the 
validity of which need testing on a case-by-case basis (see 
below). 

The diffraction contribution to the Green's function (red 
dashed line in Fig. 15) is 



PSCA/-./ 



p ikL(f — i37r/4 

(r',r,fc) = — „(()' ft k d) y 



'2nkL(f\b) 

e ikL(b,r)-i3ir/4 



G 



'2irkL{b,r) 

= G SCA (f',6,fcM<9',6>,M) x 

G SCA (b,r,k). (A3) 

We refer to this term as the (first-order) pseudo-path 
scmiclassical contribution, because it corresponds to a 
classically forbidden path. We note the different scal- 
ing of G PSCA and G rcfl w ith k. The amplitude of G rofl 
scales as 1/ y/k{L\ + L-i) while the diffractive contribu- 
tion scales as \j\/WTj\L^. At large distances (far-field) 
and/or large k, the geometric reflection amplitude dom- 
inates over the diffraction amplitude, as expected. 
The diffraction coefficient v(0' ,9,k, d) as a function of the 
incoming and outgoing angles has been calculated in the 
past in the Kirchhoff (KDA) and Fraunhofer diffraction 
approximations (FDA). Both approaches have in com- 
mon that the amplitude of the cylindrical diffractive wave 
emanating from the orifice is determined from an inte- 
gration over the orifice using as a source the amplitude 
and the phase of the unperturbed incoming wave in the 
absence of the boundary. In other words, they are im- 
plementations of Huy gens' principle according to which 
each point of the lead opening is the source of an out- 
going circular wave. For details, we refer the reader to 
Ref. 10 (for the KDA) and to Ref. 12,14 (for the FDA). 
Keller's geometric theory of diffraction (GTD) 32 has a 
different point of departure: it originates from the far- 
field approximation of the exact Sommerfeld's solution 
for the wave scattering at a wedge 58 [see Fig. 16 (a)]. 
In the following we discuss the application of the GTD 
and its refinement within the framework of the uniform 
theory of diffraction UTD 30 to the problem of the diffrac- 
tive scattering at the lead mouth [Figs. 16 (b) and (c)]. 
Starting point of our determination of the diffraction co- 
efficient v{6' ,6,k,d) is the decomposition of the orifice 
into two wedges with an inner angle of 7t/2 and outer 
angle of 37r/2 [see Fig. 16 (a)]. In the far-field limit, the 
incident wave can be regarded asymptotically as a plane- 
wave. The GTD describes the scattering of a plane wave 
at an infinitely sharp wedge 32 [Fig. 16 (a)]. The total 
wavefunction is a sum of the incoming plane wave, the re- 
flected plane wave and an outgoing cylindrical diffracted 
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Eq. (A3) but with the GTD reflection coefficient 
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FIG. 16: (Color online) (a) Diffraction at a wedge, (b) 
Diffraction at a lead described as diffraction at two wedges. 
The reference path is marked by a red dashed line. The angle 
is counted positive for a path lying on the left of the verti- 
caf axfs, for a path on the right the angle Is negative. (In 
the present example 6 is negative and 0' is positive). The 
points f and r ' are assumed to be in the far fiefd region, (c) 
Higher-order corrections enter via paths scattered between 
the wedges. 



wave, emanating from the edge. There are two discon- 
tinuities: at the shadow boundary [</>' — <f> = tt, for the 
definition of the angles see Fig. 16 (a)] and the bound- 
ary of geometric reflection ((f)' + <j> = 7r). Outside a small 
region around these two angles, the diffracted wave can 
be described by an outgoing cylindrical wave modulated 
by a smooth diffraction coefficient. The diffractive part 
of the Green's function between the two points r and r ' 
in Fig. 16 (a) can thus be approximated as: 65 



G PSCA (r', 



with 



,k) 



G SCA (r',f ,k)^D( 



SCA 



ro,r,k), 
(A4) 



D( 



' TV 



JV 



N 



+ 
N 



(A5) 



where TV = 3/2 is the exterior angle (in units of tt) of a 
perpendicular wedge and ro is the position of the corner 
of the wedge [Fig. 16 (a)]. 

We now consider the lead opening as being composed of 
two wedges [Fig. 16 (b)]. The obvious conceptual diffi- 
culty iies in the fact that the two wedges are, in generai, 
not in the far-field limit (kd ^> 1) of each other. With 
this caveat in mind, diffraction at the lead can be con- 
sidered within the GTD to result from the interference 
of two paths that are diffracted at the corners of the 
two wedges limiting the lead. Summing up the diffrac- 
tion weights and phases from the two paths, we obtain 
the Green's functions of the reference pseudo-path as in 



,GTD inl 



1 



6,k,d) = -D L {9',9)e 



I n\ — ifc^(sinfl'-l-sine) 



-D R (9',9)e +lk 2 < - sine ' +s[nf>) . 



(A6) 



In Eq. (A6) we have neglected the difference of the incom- 
ing angles at the left and at the right wedge, respectively, 
i.e., we set 6± = 62 = [Fig. 16 (c)]. Likewise, for the out- 
going angles, we set 9[ = 9' 2 = 0' . The phase differences 
of the right /left pseudo-path with respect to the reference 
path emanating from the center of the orifice can then be 
written in linear approximation with respect to the trans- 
verse iead coordinate as kAL = ±ik £ (sin 6 ; + sin 6). The 
coefficients in Eq. (A6) are defined as 



D L (9\9) 
D R {9',9) 



D(tt/2 
L>(tt/2 



0',7r/2 
0',7r/2 



9). 



(A7) 



We note that both Dl and Dr have a singularity at the 
reflection boundary for 9' — —9. However, in the sum 
of the two terms, the singularities cancel out and the re- 
sulting backscattering amplitude is perfectly smooth (see 
Fig. 17, red chain dotted line). 

The comparison of the GTD with the exact quantum 
calculation for the diffraction coefficient, \v(9, 9', k, d)\ 2 
(Fig. 17), reveals sizeable deviations. Two deficiencies 
are noteworthy: the almost complete missing of the 
back-reflection peak and the failure at grazing angles 
6' ,9 -^y^ tt/2. The diffraction coefficient should approach 
zero in this limit but, instead, converges toward a finite 
value (Fig. 17). This behavior originates from treating 
the diffraction at the iead as two independent local phe- 
nomena of diffraction at two separate wedges. 
The diffraction at a lead can be treated within the 
UTD by using a double-wedge diffraction coefficient (see 
Ref. 59,60 and a more recent paper for arbitrary configu- 
rations of the wedges, Ref. 61). The double- wedge diffrac- 
tion coefficient cannot be separated into a sequence of 
single-wedge diffraction coefficients and contains rather 
involved mathematical expressions such that the beauty 
and structural simplicity of a semiclassical approach is 
lost. We present in the following an ansatz for double- 
wedge diffraction which assumes a separation of the 
diffraction process into a sequence of diffraction events 
and we verify its vaiidity by comparing with quantum 
mechanical results from Ref. 10. We show that the short- 
comings of Eq. A6 can be (to a large extent) remedied by 
taking diffractive paths of higher order into account, i.e. 
paths that pass between the two wedges once or several 
times. This drasticaily improves the agreement with the 
quantum mechanical result. (We point to the conceptual 
similarity of our approach to the treatment of double- 
wedge diffraction illuminated by transition region fields 
by a sum over higher-order diffracted fields. 62 ) In a first 
step, we include paths that scatter once between the two 
edges [see green line Fig. 16 (c)]. There are two such 
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FIG. 17: (Color online) Absolute square of the diffraction 



coefficient \v GTD (8', 9, k, d)\ 2 (or v 



GTD-UTD 



compared with \v 



GTD-UTD 



(9',d,k,d)\ 2 for j n 



with Jn 

1, jn 



0) 

= 5 



and jmax = 10 (for the definition of j ma x see text), k = 2.5tv/cI. 
The quantum mechanical (QM) result is taken from Ref. 10. 



paths. One approaches the right wedge at angle 9. It 
is scattered into the angle ir/2 (with respect to the sur- 
face normal), and at the left wedge, it is scattered into 
the angle 9'. The other path is scattered from the left 
wedge to the right one with the same entrance and exit 
angles. The weights of this pair of paths cannot be deter- 
mined from the GTD diffraction coefficients [Eq. (A5)]: 
The GTD diffraction coefficient fails in the limit of <f> — > 
and (/>' — ► 7r [definition of angles as in Fig. 16 (a)], because 
this is in proximity of the shadow boundary, into which 
the horizontal paths are scattered. This problem can be 
overcome by invoking the uniform theory of diffraction 
(UTD). 

Contrary to the GTD, the UTD is also valid on the zone 
boundaries. The outgoing cylindrical wave is multiplied 
by a diffraction coefficient which depends not only on the 
two angles <// and <f> but also on the distances r' and r, 
and on the wavenumber k, 



UTD / if 



D 

u cot 



/ e T 

<t>,r ,r,k) = — — x 

7r + 77(0' — crq 
2N 



F ( fcJZLo^' - <r<f>)), (A8) 
\ r + r J 

where a±((3) = 2 cos 2 ^ 27rN ™ an( f n ± [ s the integer 

which most closely satisfies 2irNn ± — j3 = ±tt. The func- 
tion F is defined as a generalized Fresnel integral: 



F(x) 



00 

dre 1 



and has the asymptotic form 



F(x) = 1 



3 1 



1 



2x 4 x 2 



(A9) 



(A10) 



For x — > 00, i.e., for large distances and outside the tran- 
sition zones, F(x) = 1, and D UTD (</>', 0, r', r, k) reduces 
to D GTU ((j)' , (f>). Within the transition zone, the distance 
dependence of D UTD leads to a deviation of the scattered 
wave from a purely cylindrical wave. This is necessary 
to ensure the continuity of the total wave-function at 
the zone boundaries. Furthermore, by using the UTD 
diffraction coefficient we partially take into account the 
fact that the two wedges are not in the far-field region 
with respect to each other. 

Since the GTD fails for large (near grazing) angles, we 
opt for a piecewise construction: We combine GTD and 
UTD referred to in the following as GTD-UTD such 
that the diffraction of the path with the smaller (abso- 
lute value) of the angle is treated by the GTD and the 
path with the larger angle on the level of the UTD. Ac- 
cordingly, to first order the GTD-UTD correction to the 
diffraction coefficient can be written as 

w GTD-UTD,l ( ^ M) = 
1 p ikd(— sin #'+sinfl)/2 i 

-U L (0', -tt/2, d, k) == M+n/2, 9) 

* VZTTkd * 

1 p ikd{-\- sin 8 — sin#)/2 i 

+ -U R (6',+ir/2,d,k) == -D L (-ir/2,9), 

£ V2irkd £ 

(All) 

where [in analogy to Eq. (A7)] we have defined the UTD 
diffraction coefficients at the left (L) and right (R) wedge 
as 



U L {9',9,r,k) 

u R (e',e,r,k) 



D 
D 



UTD 
UTD 



(tt/2 - 9', tt/2 



oo, r, k) 



{tt/2 + 8', tt/2 + 9,r' -> oo, r, k). 

(A12) 



In Eq. (All), the diffraction of the incoming path (with 
angle 9), is treated on the GTD level and the diffrac- 
tion of the outgoing path (angle 6') on the UTD level. 
This construction introduces a first-order discontinuity 
("kink") at \9\ = \9'\. This kink is, however, negligible 
for small angles and visible only at large angles close to 
7r/2, in other words the formula breaks down in the limit 
6,6' — > tt/2, sec, e.g., Fig. 5 (c). (One could alterna- 
tively ignore the fact that the UTD is not multiplicative 
in the near-field and employ for both scattering events 
UTD, thus avoiding the kink. This ansatz, however, 
breaks down similarly for grazing incidence: the diffrac- 
tion coefficient docs not approach zero for the outgoing 
angle 6' — ► 7r/2. For small and medium angles this ap- 
proach behaves equally well as the presented GTD-UTD 
approach.) At large incident and/or outgoing angles the 
diffraction coefficient is already strongly suppressed such 
that the resulting error is small. 

The inclusion of the first order GTD-UTD correction 
(jmax = 1) see Fig. 17) already considerably improves the 
agreement with the quantum diffraction pattern. Higher 
order diffraction corrections include paths that are scat- 
tered several times between the wedges. This includes 



18 



paths that are incident and backscattered at an angle 
±7r/2 at the wedge. This is exactly on the reflection 
boundary. In this limit, Rcichc has shown 63 that the 
diffraction pattern of a plane wave with unit amplitude 
incident on the wedge with an angle of ±7r/2 reduces to 
a reflected plane wave with amplitude 1/2 and a cylin- 
drical wave. The Green's function of a higher order path 



which scatters j times between the wedges is, therefore, 
a product of the GTD diffraction coefficient, the UTD 
diffraction coefficient, and the Green's function for free 
propagation along the distance jd, acquiring a factor 1/2 
and a phase of tt for each reflection at a wedge. Summing 
the diffraction corrections up to order j max we obtain 



J 



.GTD-UTD ini 



,k,d) = v^ iU [e',9,k,d) + 



- ^(^^-^/2,Jd,%,(fc)e ^ *-+^(+V2,e) + ^(e^+7r/2,Jd,%,(fc)e ^ *+-^ L (-V2, 



odd: j=l 



J2 %(e^+V2,jd,%,(fc)e ^ * ++ ^i ^ (W2,0) + ^(0^-V2,.M% J (^e ^ *--^L(-7r/2,l 



even: j — 1 



(A13) 



where 



1 1 



J(kjd+(j-l)Tr) 



y/2irkjd23- 



(A14) 



and 



ef / x«/ 




k-{± sin (9' ± sin < 



(A15) 



FIG. 18: (Color online) A lead of width d coupled to a half- 
infinite plane. An incoming wave in mode m can be separated 
into two rays 9 — ±0 m = ± arcsin rrnv/dk which diffractively 
scatter at the lead wedges. The reference path is denoted by 
a red dashed line. 



Fig. 17 demonstrates that the diffraction coefficient 
v(9',0, fc,d) GTD_UTD is converged for j max = 5. Fur- 
thermore, the condition v(9' ,9,k,d) — > for 9,9' — ► 7r/2 
is fulfilled and the agreement with the fully quantum 
mechanical backscattering weight is excellent. The 
agreement deteriorates for increasing entrance angles 
9 but is still satisfying compared to the simple Fraun- 
hofer diffraction approximation (FDA), (sec Fig. 5). 
Fortunately, large angles do not play an important role 
because the overall diffraction weight is very low. The 
GTD-UTD provides a remarkable compromise between 
simplicity and accurate representation of quantum 
mechanical results for diffraction at a lead attached to a 
semi-infinite half-plane. 



2. The GTD-UTD for coupling from lead modes 
into the cavity 

The flux-normalized wave-function for mode m in a 
lead of width d oriented parallel to the x-axis is 



i>m{x,y) = 



sin (k v , m y) 



i(kx, m x-ky, m yh 



(A16) 



where k y ^ rn = mir/d is the transverse momentum compo- 
nent and k x . m = yfk 2 — (mir/d) 2 the longitudinal com- 
ponent. Eq. (A16) can be viewed as two rays emanating 
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FIG. 19: (Color online) The absolute square of the coupling coefficient \c m {6,k,d)\ 2 multiplied by k Xtm for m — 2 at two 
different wavenumbers (a) k = 2.1n/d and (b) k — 2.5ir/d within the GTD and the GTD-UTD. Note the breakdown of the 
GTD for 6 -» ±7r/2 near k « 2%/d. 



with angles ±6 m = ± arcsin(m7r/<ifc) (Fig. 18). In the 
GTD approximation, each ray hits an edge of the lead 
mouth and the two cylindrical waves emanating from the 
edges cause, in turn, an interference pattern at large dis- 
tances. The GTD diffraction for scattering from the lead 
mode m into the half-space (in our application the cav- 
ity) can be written in direct analogy with Eqs. A6 and 
A7 as 



The two paths have a phase difference of mir at the lead 
mouth and thus a phase difference ±rmr/2 relative to the 
reference path that starts at the center of the orifice (see 
Fig. 18). With the flux normalization factor yJk Xytn from 



Eq. (2) ^/k~c. 



GTD; 



, k, d) is dimcnsionlcss 



We include now higher-order scattering events on the 
level of the UTD in order to improve Eq. (A17), 



„GTD 



(e,k,d) 



(A17) 



with the scattering coefficients at left and right wedge 



D R (6,6 m ) = D{\ 



3tt 



(A18) 
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where 



±± 



±k— sin 6? ± V^L 
2 2 ' 



gj(k) is given in Eq. (A14) and Ul, U r are given by 
Eq. (A12). Fig. 19 illustrates that the UTD corrections 



(A20) 
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become most important when the value of k is close to a threshold, i.e., when the angle 9 m is close to tt/2. 
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